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ABSTRACT 

We develop a method for estimating the host galaxy dust extinction for type la 



supernovae based on an observational coincidence first noted by [Lira (1995)| , who 
found that the B — V evolution during the period from 30-90 days after V maximum 
is remarkably similar for all events, regardless of light curve shape. This fact is used to 
calibrate the dependence of the B max — V max and V max — I ma x colors on the light curve 
decline rate parameter Amis(-B), which can, in turn, be used to separately estimate 
the host galaxy extinction. Using these methods to eliminate the effects of reddening, 
we reexamine the functional form of the decline rate versus luminosity relationship and 
provide an updated estimate of the Hubble constant of H a = 63.3 ± 2.2(internal) =t 
3.5(external) km s _1 Mpc -1 . 
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Subject headings: distance scale - supernovae: general 



Introduction 



In recent years, considerable attention has been focused on the utility of type la supernovae 
(SNe la) as cosmological standard candles (e.g., Hamuy et al. (1996b) , Perlmutter et al. 



|(1995)| , |Ricss et al. (1996a)| , Sandage et al. (1996)| ). While there is now abundant evidence for the 
existence of a significant dispersion in the peak luminosities of these events at optical wavelengths, 
the absolute magnitudes fortuitously appear to be closely correlated with the decay time of the 
light curve ( [Phillips (1993)1 , |Hamuy et al. (1995)| , |Riess et al. (1995)| , |Hamuy et al. (1996a)| , |Riess 
et al. (1996a)|). After correction for this effect, dispersions of ~0. 15-0. 20 mag were obtained by 



Hamuy et al. (1996b) in the BVI Hubble diagrams of 26 SNe la in the redshift range 0.01 < z < 
0.1, and a Hubble constant of 63.1 ± 3.4 (internal) ± 2.9 (external) km s _1 Mpc -1 was derived 



via reference to four nearby SNe la with Cepheid-calibrated distances. Although Hamuy et al 



(1996b) applied a color cut to their sample to eliminate the reddest events, no corrections were 



applied for possible host galaxy dust extinction of the supernovae in either the distant or nearby 
samples. Not only might this affect the derived slope of the decline rate versus luminosity relation, 
but the correctness of the Hubble constant depends on the assumption that the mean reddening 
of the 26 distant SNe la in the Hubble flow is essentially the same as that of the four nearby 
calibrating SNe. 



In a parallel study which included many of the objects observed by |Hamuy et al. (1996b) , 



Riess et al. (1996a) developed an empirical method which they called "Multicolor Light Curve 



Shapes" (MLCS) to estimate simultaneously the luminosity, distance, and total extinction of 
each SN la. MLCS employs linear estimation algorithms to create a family of SN la light 
curves and color curves, and offers considerable elegance in error modeling. However, this first 



implementation by Riess et al. (1996a) was based on essentially the same "training set" of only 



nine nearby supernovae used by Phillips (1993) which suffered from uncertain extinction estimates 
and inhomogeneous secondary distance indicators. This training set was sufficient to determine 
the existence of the relations between SN la light curve shape, color, and luminosity, but not ideal 
for determining the precise values of these relations. Indeed, a recent reanalysis of MLCS in the 



B and V bands by [Riess et al. (1998)| using a much larger training set as well as redshifts as 
the distance indicator to determine the SN la luminosities improves the precision of the method. 
Additionally, this reanalysis reduces the strong relation between light curve shape and B — V color 
at maximum light initially found by Phillips (1993) and Riess et al. (1996a)| . 



In this paper, we seek to establish an independent procedure for determining host galaxy 



reddenings for SNe la using an observational coincidence first noted by Lira (1995) . This method 
is developed in § 2 and used to estimate the intrinsic B max — V max and V max — I ma x colors for a 



-3 - 



sample of ~60 well-observed SNe la. We then apply this information in § 3 to derive reddening-free 
relations between the decline rate parameter Amis(B) (Phillips (1993)) and the peak absolute 
magnitudes in BVI, and to update our estimate of the value of the Hubble constant. 



2. Color Excesses and Intrinsic Colors 

2.1. E(B - V) Tail 

In comparing the color evolution of several apparently unreddened SNe la representing the 



full range of observed light curve decline rates, [Lira (1995)| found that the B — V colors at 
30-90 days past V maximum evolved in a nearly identical fashion. This property, which was 
independently discerned by Riess et al. (1996a) , is illustrated in Figure [j] where we plot the 



B — V color evolution from CTIO photometry of six SNe la covering a wide range of decline rates 
(0.87 < Ami5(.B)< 1.93) and which likely suffered little or no reddening from dust in their host 
galaxies. The likelihood of a supernova being unreddened was judged by three basic criteria: 1) 
the absence of interstellar Na I or Ca II lines in moderate-resolution, high signal-to-noise spectra, 
2) the morphology of the host galaxy (SNe that occur in E or SO galaxies are less likely to be 
significantly affected by dust), and 3) the position of the supernova in the host galaxy (SNe lying 
outside the arms and disks of spirals are more likely to have low dust reddening). As discussed by 
Lira, it is not surprising that the colors of SNe la at these late epochs are so similar since their 
spectra also display an impressive uniformity. This is precisely the epoch where the Fe-Ni-Co core 



dominates the spectrum which is rapidly evolving into the "supernebula" phase (see Phillips et al 



19921) 



From a least squares fit to the photometry of four of the SNe displayed in Figure (1992A, 
1992bc, 1992bo, and 1994D), Lira derived the following relation to describe the intrinsic B-V color 
evolution: 



(B-V) = 0.725 - 0.0118(£y -60), (1) 

where ty is the phase measured in days since V maximum. This fit, which is plotted in Figure 
is valid over the phase interval 30 < ty < 90, and can be applied to any SN la with BV coverage 
extending to at least ty = 30 days to derive an estimate of the color excess E(B — V). Note in 
Figure |l| that individual SNe can display systematic residuals with respect to the Lira fit (i.e., the 
points lie mostly bluer or redder than the line) suggesting that there is an intrinsic dispersion 
about equation |l|. This dispersion amounts to 0.04 mag for the four SNe used by Lira, and 
0.06 mag for the six SNe shown in Figure H; in calculating color excesses from equation |], we shall 
adopt a value of 0.05 mag in our error analysis. 

For relatively bright SNe observed with CCDs, the precision of the photometry obtained at 
~l-3 months after maximum is typically still quite good (0.02-0.05 mag). However, for more 
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distant events, the late-time photometry coverage tends to be relatively sparse with typical errors 
~0.1-0.2 mag for individual observations, and values of E(B — V) derived directly from equation [j] 
will often have large associated errors. For such SNe, a better technique is to use template fits of 



the type employed by Hamuy et al. (1996c) to derive an estimate of the observed B — V color 
at some fiducial epoch in the interval 30 < ty < 90. This value can then be compared with 
the intrinsic color predicted by equation [l] for the same epoch to derive the color excess. This 
procedure has the advantage of using the entire set of BV photometry to estimate a late-epoch 
color rather than just the subset of data obtained at 30 < ty < 90. 

For a sample of 62 well-observed SNe la with z < 0.1 consisting of 1) the "Calan/Tololo" 
sample of 29 SNe la published by Hamuy et al. (1996c)| , 2) 20 SNe la published in |Riess et al] 



(1999) which we shall refer to as the "CFA" sample, and 3) 13 nearby, well-observed SNe la 
(1937C, 1972E, 1980N, 1981B, 1986G, 1989B, 1990N, 1991T, 1991bg, 1992A, 1994D, 1996X, 
and 1996bu — see Table |] for photometry references), we have calculated color excesses from 
equation |l| via the two techniques described above — i.e., both directly from the late-time 
photometry and via light curve template fits. In the following, we shall refer to these as the 
"Direct" and "Template" methods, respectively. For the "Template" method, we chose ty = 53 
days as the fiducial reference epoch for which equation [l] predicts (B — V) = 0.808. To calculate 
the host galaxy reddening, E(B — V)f los t, the observed supernova color evolution must be corrected 
for both Galactic reddening and the effect of redshift (the "K correction") via the equation 



(B - V) c 



(B - V) obs - E(B - V) Gal - Kb-v, 



(2) 



where (B — V) Q b s is the observed color, E(B — V)c a i is the Galactic reddening, for which we use 
the values from Schlegel et al. (1998) , and Kb-v is the K correction. The host galaxy reddening 
is then given by 



E(B - V) host — (B — V) c 



(B-V) . 



(3) 



As discussed by Rless et al. (1998") and Nugent et al. (1999) , Kb-v is not only a function 
of redshift but also the SN color. Since we are concentrating on an epoch where all SNe la 
show essentially the same B — V color, we need only worry about color differences arising from 
different amounts of host galaxy reddening. In principle, an iterative procedure should be followed 
where Kb-v for zero host galaxy reddening ( Hamuy et al. (1993) ) is taken as a first guess, 
E{B — V)host is calculated via equations [2] and |3|, and the Kb~v correction is then adjusted if 
E(B — V)h os t is found to be significantly greater than zero. However, for the range of host galaxy 
reddenings (E(B — V) < 0.8) typical of our sample, the variation of Kb~v due to this effect is 
only 0.00-0.02 mag. This is much smaller than the precision of our reddening indicator, and so we 
therefore adopt the K#_y corrections for zero extinction with no iteration. 



As pointed out by 



Kim (1997) 



, the observed color excess of a SN la will vary slightly with 
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light curve phase due to the significant color evolution that these events undergo. There will 
also be a generally smaller dependence on the total amount of dust extinction itself (e.g., see 
Blanco (1956)). Using the IRAF task "deredden", which employs a standard Galactic interstellar 



reddening curve flCardelli et al. (1989) ), and a program which we have developed for deriving 



synthethic photometry from spectra, we have calculated observed color excesses as a function of 
both the light curve phase and the "true" reddening, E(B — V)t rU e, for a database of more than 
30 SNe la spectra in order to characterize and correct for these effects. Note that E(B — V)t rue 
parameterizes the amount of extinction along a particular line of sight, and not the observed 
color excess which, as stated, will vary as a function of phase and reddening. To convert from 
E(B — V) \ )S to E(B — V)true at ty = 53 days, our reference epoch for the "tail" evolution phase, 
we find the following approximate relation to be valid: 

E(B - V)true = 1.018/(1/^(5 - V) obs - 0.072). (4) 

The observed decline rate of a SN la will also be a weak function of the dust extinction which 
affects the light curves. This follows from Figure H where we have used the above-described 
calculations to examine the variation of 

R B = A B /E(B - V)true (5) 

with light curve phase for three different values of E{B — V)t rue . In this equation, Ab is the 
observed extinction which is equal to the difference between the synthetic B magnitudes calculated 
for a reddened and unreddened spectrum. Combining equation [5] with the definition of the 
observed decline rate parameter, 

Ami 5 (B) obs = B obs (+15 days) - B obs (max), 

gives 

Am 15 (B) obs = Am 15 (B) true + [R B (+15 days) - R B (max)]E(B - V) true 
From Figure H this can be shown to reduce to the following approximate relation: 

Am 15 (B) true ~ Am 15 (B) obs + 0.1E(B - V) true (6) 
Note that reddening acts to decrease the decline rate, which is opposite to the conclusion of 



Leibundgut (1988) who employed a formula from gchmidt-Kaler (1982)| which gives Ab as a 



function of the intrinsic color, (B — V) a , and the observed color excess, E(B — V) a b s . Although 
this formula was derived for normal stars, it does a surprisingly reasonable job of predicting the 
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shape of the variations of As with light curve phase implied by Figure |2[ However, to use this 



formula correctly, one must take into account that E(B — V) f, s also varies with phase; Leibundgut 



assumed that E(B — V) b s was constant with phase, which led him to incorrectly conclude 
that the effect of reddening was to increase the decline rate. 

Figure [| illustrates application of the "Direct" and "Template" methods to the photometry 
of two representative SNe observed in the course of the Calan/Tololo survey. Comparison of the 
E(B — V)host values (converted to "true" reddenings via equation ||) derived via both techniques 
for 54 of the SNe in our sample is shown in Figure £§ In general, there is good agreement between 
the two reddening determinations, with no evidence for a difference in zero points (the weighted 
average of the difference amounts to 0.01 ± 0.01). The most discrepant points correspond to 
SNe 1992ag and 1996ai (not plotted). SN 1992ag, which was discovered and observed in the 
Calan/Tololo survey, occurred near the center of an edge-on spiral host galaxy, and we suspect 
that the B photometry obtained at late epochs was systematically affected by the difficulty of 
subtracting the host galaxy. A similar explanation may apply to the severely reddened SN 1996ai 
for which we find E(B — V)Tempiate = 1-97 ± 0.08 and E(B — V)Direct = 2.43 ± 0.19; in this case, 
the E(B - V) Direct estimate is based on a single color measurement which we suspect may be 
more uncertain than the formal photometry errors would suggest. 

The zero point of these color excesses can be checked by examining the reddenings obtained 
for the subset of SNe la with E/S0 host galaxies, where we expect there to be little or no host 
galaxy dustf]. For the 21 SNe la in our sample with E/S0 hosts, we defined E(B — V)to,u to De 
the weighted mean of the host galaxy color excesses E(B — V)Tempiate and E(B — V)Directf^ We 
find an average value of E(B — V)xaii = 0.05 ± 0.03, which decreases to 0.03 ± 0.03 if we limit the 
sample to the eight SNe in E galaxies. We do not consider either of these results to be significantly 
different from zero. Moreover, there are additional arguments for retaining the Lira zero point: 

1. The SNe used to determine the Lira zero point were all relatively nearby. We therefore 
have precise morphological classifications for the host galaxies. Being nearby, it's also easier to see 
signs of dust directly in the host galaxies. In addition, we have high signal-to- noise spectra which 
allow us to put limits on the amount of interstellar gas present in the line of sight. Although Lira 
used only four SNe for her zero point, the six shown in Figure [l] give the same zero point to 0.01 



3 Although there is some evidence for a diffuse component of dust in nearby E galaxies that would be difficult to 



detect directly (e.g., see Goudfrooij fc de Jong (f995) ), the SNe with E hosts in our sample occurred at relatively 



large projected distances (> 10 kpc) from the centers of their host galaxies where the maximum amount of reddening 



consistent with observed color gradients is typically E(B — V) < 0.05 ( Wise fc Silva (1996) ). Thus, it seems reasonable 



to expect that SNe in E galaxies should give a reliable determination of the reddening zero point. 

4 ln taking a weighted mean of E(B — V)oirect and E(B — V)Tem P iate, we have assumed these to be independent 
measurements of the reddening. Technically speaking, this is not fully the case since the tail photometry used to 
measure E(B — V) Direct is a subset of the entire phoometry set which we use in calculating E(B — V)Tem P iate- 
However, for a typical SN in our sample, the tail data amount to less than half of the full set of photometry and 
enter into the template fits with low weight due to their larger measurement errors. 
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mag. 

2. Of the 21 SNe in our sample with E/SO host galaxies, many are relatively distant. Some of 
the morphological classifications of the hosts could therefore be in error (e.g., at large distances 
it's sometimes difficult to distinguish an SO galaxy from an Sa). Hence, we can't be certain that 
some of these galaxies don't have significant dust. Also, very few of the spectra obtained of these 
SNe have sufficient signal-to-noise to look for interstellar absorption lines. 

Thus, we adopt the zero point of equation |l| without modification. In any case, it should 
be recalled that the precision of any single measurement of E(B — V)to,u wu 1 be limited by the 
intrinsic dispersion about equation |l| which, as stated earlier, appears to be ~0. 04-0. 06 mag. 



2.2. E(B - V) Max and E(V - I) M ax 

The "tail" evolution of the intrinsic B — V color as an indicator of host galaxy dust reddening 
offers the advantage of being insensitive to the initial decline rate (or luminosity) of the SN, but 
the obvious disadvantage of this method is that it requires the light curves to be followed to 
several magnitudes below maximum light, which in practice is often impossible for distant SNe la. 
Ideally, we would prefer to be able to use the maximum-light colors to determine the host galaxy 
reddening. As is seen in Figure || (and pointed out by Lira (1995) ) , the B — V colors at maximum 



light of most SNe la show a fairly small dispersion, although the fastest declining events (such as 
1991bg) are clearly significantly redder. In principle, we can calibrate any dependence on decline 
rate by using the "tail" method to identify a sample of unreddened SNe la representing a wide 
range of initial decline rates. 

Columns 2-4 of Table ||| list the measured decline rate parameter Amis(i?) obs , the color excess 
E(B — V)cai due to Galactic reddening, and the E(B — V)tuU values calculated as per § 2.1 for the 
62 SNe la in our full sample. We find that a total of 20 events have values of E(B — V)xaU < 0-06 

- i.e., less than or equal to the standard deviation of a single reddening measurement. In the 
following, we shall assume that these SNe la were essentially unreddened by dust in their host 
galaxies. The "colors" B max — V max and V max — I max for these 20 events are plotted versus the 
decline rate parameter Ami5(B) in Figure ||. Here we differentiate between those SNe la with 
light curve coverage beginning < 7 days after the epoch of B max and those for which the first 
observation was not obtained until > 7 days after the epoch of B max since we expect that the color 
excesses determined for the latter SNe will be less precise than those for events observed within a 
week of maximum light. For decline rates in the range 0.9 < Am\^(B) < 1.6, the variation in color 
is approximately linear in both colors. The fastest-declining events, 1991bg and 1992K, clearly 
do not follow this trend, presumably due to the appearance in their spectra of strong absorption 
features of low-ionization species such as Ti II as a consequence of lower effective temperatures 



(see [Nugent et al. (1995)| ). 



Over the linear portion of the relations shown in Figure ||, we calculate the following weighted, 
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least squares fits to the SNe la with light curve coverage beginning < 7 days after B„ 



(Br, 



V„. 



-0.070(±0.012) +0.114(±0.037)[Am 15 (5) - 1.1], a = 0.030, n = 14 (7) 



-0.323(±0.017) +0.250(±0.056)[Am 15 (B) - 1.1], a = 0.042, n = 11 (8) 



These fits are plotted as solid lines in Figure |||. For reference, we also show the relations derived by 
Phillips (19931 , |Hamuy et al. (1996a)| , the original formulation of MLCS ( [Riess et al. (1996a)[) and 
the revised version of [Riess et al. (1998)] . Clearly the slopes of the B max — 

Vmax relations found by 



Phillips (1993) and the original MLCS are too steep, whereas both the Hamuy et al. (1996a) and 



the new MLCS provide much better approximations to the data. The small offset between our fit 
and that of Riess et al. (1998) is due in part to our usage of the |Schlegel et al. (1998)| Galactic 
reddenings, which are on average ~ 0.02 mag greater than the Burstein fc Hcilcs (1982) values 
employed by Riess et al. (1998") ; the larger offset with respect to the Hamuy et al. (1996a)| fits is 
due to the same effect, and also to the fact that Hamuy et al. (1996a) did not correct their data 
for host galaxy extinction. Curiously, the Phillips (1993) relation for V max — I max is a poor fit to 
the data, whereas the original MLCS reproduces the dependence quite well. 

To use the relations given by equations § and ^ to estimate host galaxy reddenings, we 
procede in a fashion analagous to that described in § 2.1. The observed color excess is defined as 



E(B — V)host — {B m ax — V r 



max Icorr 



{.Bmax Vrr. 



(9) 



where (Bmax - V, 



max "max Icorr 



is given by 



(Br, 



Knai)t 



(Br, 



V, 



max lobs 



E(B - V) 



Gal 



K B 



max 1 rn ax 



(10) 



In principle, since Ami^(B) obs is a function of the reddening, equations ^ and [T(] should be 
employed in an iterative fashion, but the slopes of equations [/] and || are sufficiently shallow that 
this effect can be ignored for all but the most heavily reddened events. The observed color excesses 
are then converted to "true" color excesses via the following approximate relations: 



E(B - V) 



true 



0.981/(1/^(5 - V) ohs - 0.050) 



(11) 



E(V - I)true = 0.989/(l/£(F - I) obs - 0.004) 



(12) 
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Note that these color excesses, which we shall refer to as E(B — V)mo,x and E(V — I) Max, are valid 
only for SNe la with decline rates in the range 0.9 < Ami^^B) < 1.6 which also show "normal" 
spectra at maximum light - i.e., with strong Si II A6355 absorption. In particular, this method 
should not be used to calculate reddenings for SN 1991T-like events which show very different 



spectra at maximum and therefore may have dissimilar intrinsic colors (e.g., see Phillips et al 



i L992) 



Values of E(B — V) Max and E(V — I) Max for the SNe in our sample are listed in columns 
5-6 of Table [|, and are plotted versus E(B — V)i- a u in Figure ||. Except for SN 1992ag, which 
we have already discussed, the agreement between E(B — V)mux and E(B — V)to,U 1S generally 
excellent. The relationship between E(V — I) Max and E{B — V)TaU, while clearly significant, 
is not as impressive. Perhaps this is not surprising in view of the larger dispersion observed in 
the E(V — I)mox versus Ami^(B) relation for relatively unreddened SNe la (see Figure || and 
equation ^) which we suspect arises from two main causes: the generally poorer sampling of the 
/-band light curves for many of the Calan/Tololo SNe, and the more complicated nature of the 
shapes of the I-band light curves which makes predicting I max from template fitting challenging. 
Nevertheless, we consider the V max — I m ax color to be a useful alternative method of identifying 
significantly reddened SNe Ia0. 

In column 7 of Table |2] we list the average color excess E(B — V)Avg which we take to be the 
weighted mean of E(B — V)Taih E{B — V)mox, and 0.8E(V — I) Max- The multiplicative factor of 
0.8 for E(V — I) Max corresponds to an assumption of a standard Galactic interstellar reddening 



curve (|Dean et al. (1978) ; see also § 3). In addition, to handle negative reddening values we have 



applied the "Bayesian filter" used by Riess et al. (1998) which assumes a one-sided Gaussian a 



priori distribution of As values with a maximum at zero and a = 0.3 mag. For the subsample of 
19 SNe la in E and SO host galaxies with 0.9 < Amis (5) < 1-6, we find 



E(B - V) A vg = 0.02 ± 0.03, a = 0.05, n = 21, 

which implies that the zero points of equations [l], and |8] are reasonable, and that the typical 
precision of a measurement of E(B — V)Avg for any single SN la is ~0.05 mag. As expected, 
we find that the most reddened SNe were observed in spiral galaxies: 82% of the SNe with host 
galaxy reddenings greater than 0.05 mag occurred in spirals, whereas only 29% of the SNe with 
E(B — V)Avg < 0.05 mag had spiral galaxy hosts. 

Figure ^ shows the total reddening, E(B — V)avqi derived for the Calan/Tololo SNe la plotted 
as a function of host galaxy redshift. From the Hubble diagram fits given in § 3 and the definition 



5 Strictly speaking, color excesses estimated via equations ^ and ^ are not fully independent of the E(B — V)T a u 
measurements for the same SNe since a change in the zero point or slope of the Lira relation (equation |l|) could 
change somewhat the exact sample of unreddened events used to derive equations |^ and |^. Also, an uncertainty in 
Vmax will produce a correlated error in E(B — V)uax and E(V — I) Max since both use the same measurement. 
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of the observed magnitude of a supernova as B b s = B Q + As, it can easily be shown that the 
maximum allowed color excess for a SN as a function of redshift in a magnitude-limited search is 

E(B -V)~ (B UmU - Uog{cz) + 3.57)/4.1. (13) 
The dashed line in Figure ^ shows equation 13 plotted for the limiting magnitude of the 



Calan/Tololo search, Bn mit ~ 19 ( Hamuy et al. (1993) ). This curve provides a good fit to the 



upper envelope of the observed reddenings, particularly at higher redshifts. The fact that SNe in a 
magnitude-limited search will have preferentially smaller dust reddenings as a function of redshift 
will naturally lead to a greater percentage of SNe being discovered in E/SO galaxies at higher 
redshifts, exactly as is observed in the Calan/Tololo sample ( [Hamuy et al. (1996b)| ). These results 
and predictions provide confidence in the validity of our techniques for measuring the reddenings 
of SNe la. 



3. The Reddening-Free Decline Rate versus Luminosity Relation 

Once the reddening of a SN la is known, correcting the observed maximum-light magnitudes 
for this extinction is straight-forward. The general formula for a given filter is 

m = (m obs - A m (Gal)) - K m - A m (host), 

where m & s is the observed maximum- light magnitude, A m {Gal) is the Galactic component of 
extinction, K m is the K correction (which, as noted in § 2.1, is a function of both redshift and 
reddening), and A m (host) is the host galaxy component of extinction. From the data shown in 
Figure |2] and similar calculations for the V and / bands, the dependence of the extinction on color 
excess at maximum light can be approximated as 

A B = [4.16 - 0ME(B - V)true]E{B - V) trU e, (14) 
Ay = [3.14 - 0.02E(B - V) true ]E(B - V) true , (15) 

and 



Equations 
extinction. 



14 16 



Ay = [1.82 + 0.01^(5 - V)true]E(B - V) true . (16) 

are valid for calculating either the Galactic or host galaxy components of the 
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The left half of Figure || shows the absolute magnitudes corrected only for Galactic reddening 
for the full sample of 29 Calan/Tololo and 12 CfA SNe with z > 0.01 plotted as a function of the 
observed decline rate parameter Am 15 (i?) ofes . The distance to each SN was calculated from the 
redshift of the host galaxy (in the CMB frame) and an assumed value of the Hubble constant of 
H = 65 km s _1 Mpc _1 . Following Hamuy et al. (1996a), a peculiar velocity term of 600 km s _1 



has been included in the error bars of the absolute magnitudes. 

In the right half of Figure || this plot is repeated for the subsample of 18 of these 41 SNe for 
which we find insignificant host galaxy reddening (E(B — V)Avg < 0.05). This diagram reveals 
more clearly than ever the true nature of the peak luminosity-decline rate relations for SNe la. 
Elimination of the events with significant host galaxy reddening not only decreases the scatter in 
the relations, but also shows these to be nonlinear in the B and V bands, and probably also in /. 
(Note that Riess et al. (1998)] also found this dependence to be nonlinear and used a quadratic 



term in the revised MLCS.) In the right half of Figure pL we plot for comparison as dashed lines 



the linear fits given by Hamuy et al. (1996a) for the 26 Calan/Tololo SNe with B max — V max < 0.2. 
These show best agreement in / where reddening effects are smaller, but are progressively worse 
approximations to the actual peak luminosity-decline rate relations in the V and B bands due to 
the combined effects of uncorrected host galaxy extinction and the curvature which is apparent in 
the true relations. 

For this "low host galaxy reddening" subsample of the Calan/Tololo and CfA SNe with 
decline rates in the range 0.85 < Ami^(B) < 1.70, we calculated a least square fit to a quadratic 
peak luminosity-decline rate relation for the V band data. This fit was then combined with 
equations [?] and || to give corresponding relations for the B and I bands. Table |3| lists the resulting 
parameters for these fits expressed in terms of the magnitude difference, AM max , with respect to 
the maximum light magnitude of a SN with Ami^(B) = 1.1. Note that the low values of xt are 
probably due to an overestimate of the errors by |Hamuy et al. (1996a) which we have adopted 



verbatim. The impressively low dispersions (0.09-0.13 mag) of these fits strongly support the 
validity of the methods described in this paper for estimating the host galaxy dust reddening. For 
reference, Hamuy et al. (1996a)| , who employed a color cut of B max — V max < 0.2 and assumed 



linear relations for Amis(i?) versus M max , obtained dispersions of 0.17, 0.14, and 0.13 mag for 
the fits in B, V, and /, respectively. The wavelength dependence of these values is due to a 
combination of uncorrected host galaxy extinction and the non- linear shape of the Amis (B) versus 
M max relation in B and V; when these effects are taken into account the dispersion is uniformly 
low in BVI. 

From the average values of the host galaxy reddening given in Table [2] and the decline rate 
versus luminosity relations given in Table ||, we have reexamined the corrected Hubble diagrams 
in BVI for the 26 Calan/Tololo with B max — V max < 0.2, which Hamuy et al. (1996a, 1996b) 
referred to as their "low-extinction" sample. To provide an absolute calibration for these diagrams 
and, hence, an estimate of the Hubble constant, we use the same four SNe (1937C, 1972E, 1981B, 



and 1990N) with Cepheid distances employed by Hamuy et al. (1996b) . Three different cases were 
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considered: 1) Correcting for Galactic reddening only, 2) correcting for Galactic reddening and the 
Ami^(B) versus M max relation, and 3) correcting for Galactic reddening, host galaxy reddening, 
and the Amis (5) versus M max relation. 

Averaging the results for BVI, we find the Hubble relations illustrated in Figure |9[ (Note 
that the errors given in this figure for H are internal only, and do not include the external error 
due to uncertainties in the zero point of the Cepheid period-luminosity calibration, which we 
estimate to be an additional 3-4 km s _1 Mpc -1 .) Not surprisingly, for the first two cases we 



find Hubble constants and dispersions which are similar to those given by Hamuy et al. (1996b) 



for comparable corrections. Note that case 1 is essentially the method employed by jSandage~et 



al. (1996) , who ignored the existence of the Amis(i?) versus M max relation; the lower Hubble 
constant results from the fact that the average decline rate of the 4 calibrators with Cepheid 
distances (Amis(S) = 0.97) is significantly lower than the mean value for the 26 Calan/Tololo 
SNe (Ami5(i?) = 1.24), meaning that the calibrators are on average ~ 0.2 mag more luminous 
than the Calan/Tololo SNe. Correcting for the Ami^(B) versus M max relation (case 2) eliminates 
this bias and significantly decreases the dispersion in the Hubble diagram. Applying the host 
galaxy reddening corrections (case 3) decreases the dispersion still further, illustrating the merit 
of our approach to calculating these corrections, although we find essentially the same value of 
H as when this correction is ignored (case 2). This is ascribed to the coincidence that the median 
reddening {E{B — V)host = 0.05 mag) of the 26 distant SNe happens to be virtually identical to 
the median reddening {E{B — V)host = 0.06 mag) for the four calibrating events. 

From fits to the fully corrected (i.e., for Galactic + host galaxy reddening, and the &m.\§(B) 
versus M max relations) magnitudes of our complete sample of 28 Calan/Tololo and 12 Cfa SNe 
with z > 0.01 and Amis(i?) < 1.70, we find the following relations for the Hubble constant in 
BVI: 



logH (B) = 0.2{M% ax - 0.786(±0.398)[Ami 5 (5) - 1.1] - 

0.633(±0.742)[Ami 5 (5) - l.l] 2 + 28.671(±0.043)}, (17) 

logH (V) = 0.2{Ml, - 0.672(±0.396)[Am 15 ( J B) - 1.1] - 

0.633(±0.742)[Ami 5 (B) - l.l] 2 + 28.615(±0.037)}, (18) 

logF (I) = 0.2{M^ - 0.422(±0.400)[Am 15 (£) - 1.1] - 

0.633(±0.742)[Am 15 ( J B) - l.l] 2 + 28.236(±0.035)}. (19) 



Combining equations |i~7|-|i~9| with the corrected absolute magnitudes of the six well-observed SNe la 
(1937C, 1972E, 1981B, 1989B, 1990N, 1998bu — see |Saha et al. (1999)| and |Suntzeff et al. (1999) 
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for details and references) with HST Cepheid distances gives a value of the Hubble constant of 63.3 
db 2.2(internal) ± 3.5(external) km s _1 Mpc -1 . Until more Cepheid distances become available for 
nearby, well-observed events, this value represents our best estimate for the Hubble constant. 

Finally, we briefly address the question of whether the extinction law of the host galaxy 
dust affecting SNe la is consistent with that typically observed for dust in the disk of our own 



Galaxy. Using their original formulation of MLCS, Riess et al. (1996b) concluded that this was 
the case, but since this same version of MLCS incorrectly predicted a strong relation between 
light curve shape and the B — V color at maximum light, it is worth reexamining the issue. We 
do this by plotting the absolute magnitudes of the SNe in our sample versus the B max — V max 
color. However, since both of these quantities are a function of the decline rate we must first 
correct them to a "standard" decline rate (we choose Amis (5)= 1-1) using the coefficients of the 
decline-rate-dependent portions of the relations given in Table || and equation |?]. Specifically, we 
calculate, 

M{B) U1 = M(B) max + 0.786 [Amis (B) - 1.1] + 0.633[Ami 5 (B) - l.l] 2 , (20) 
M(y) 1A = M(V) max + 0.672[Ami 5 (S) - 1.1] + 0.633 [Amis (B) - l.l] 2 , (21) 
M(I) 1A = M(I) max + 0.422[Am 15 (£) - 1.1] + 0.633[Ami 5 (£) - l.l] 2 , (22) 

and 

(B max - V m ax)l.l = {Bmax ~ V m ax) ~ 0.114[Ami 5 (B) - 1.1]. (23) 

These quantities are plotted in Figure |l^ for the subsample of Calan/Tololo and CFA SNe la with 
z > 0.01 for which distances were calculated assuming _ff =63.3, and also for the SNe la with 
Cepheid-based host galaxy distances. In both cases, the samples were restricted to those events 
with light curve coverage beginning < 5 days after B max and for which we have derived non-zero 
values of E{B — V)^ vg . Plotted in Figure |l(] for reference are the canonical Galactic reddening 
vectors with slopes of 4.1 in B, 3.1 in V, and 1.8 in /. A linear least squares fit to the data gives 
slopes of 3.5 ± 0.4, 2.6 ± 0.4, and 1.2 ± 0.4, which agree to within ~ 1 — 1.5 a with the Galactic 



values, and which are consistent with typical values of R observed in nearby galaxies (e.g., Hodge 



fc Kennicutt (1982)| , [Bouchet et al. (1985)1 , jye fc Richter (1985)1 , Prosch fc Loinger ( 1988)|) 



Hence, we conclude that the host galaxy dust is similar in its properties, on average, to the dust 
in our own Galaxy. 
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4. Conclusions 

In this paper we have presented methods for estimating the host galaxy dust extinction of 
SNe la at late times (30 < ty < 90), when the B — V evolution is remarkably similar for all 
events, and at maximum light where the colors B max — V max and V max — I m ax for normal-spectra 
SNe la display only a mild dependence on light curve shape. From this, we are able to deduce 
the reddening-free relations in BVI for the dependence of absolute magnitude on the decline rate 
parameter Amis Fits to these data yield remarkably low dispersions (0.09-0.13 mag) in all 
three colors, and lead to a revised value of the Hubble constant of H a = 63.3 ± 2.2(internal) ± 
3.5(external) km s _1 Mpc -1 . 

We thank Adam Riess for his valuable comments on this paper, and for providing the new 
MLCS fit in Figure |5[ J.M. acknowledges support from "Catedra Presidencial de Ciencias 1995" 
and FONDECYT grant No. 1980172. 
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Fig. 1. — B — V color evolution for six SNe la which likely suffered little or no reddening from dust 
in their host galaxies. These six events, whose Ami5(i?) parameters are indicated in parentheses, 
cover a wide range of initial decline rates and peak luminosities. The solid line corresponds to the 
Lira (1995) fit (equation |l]) to the color evolution during the phase interval 30 < ty < 90. The 



epoch of B max is assumed to occur 2 days before the epoch of V m ax ( Leibundgut (1988j| ) 



Fig. 2. — Variation of Rb = Ab/E(B — V)tme as a function of phase for different values of 
E(B — V)true- Note that Ab is the observed extinction in B, which is a function of the light curve 
phase. 

Fig. 3. — Illustration of the the "Direct" (upper panels) and "Template" (lower panels) techniques 
for estimating the host galaxy reddening of SNe la from the B — V evolution during the phase 



interval 30 < ty < 90. The SN photometry is from Hamuy et al. (1996c) and has been corrected 
for both Galactic reddening and the K correction. The dashed line in each graph corresponds to 
the Lira (1995) fit (equation [l]) to the color evolution of unreddened SNe la. The solid lines in 



the upper panels are the Lira relation shifted to fit the photometry of the SNe in the interval 30 
< ty < 90; E(B — V) Direct is then defined as the difference between the solid and dashed lines. 



In the lower panels, the solid lines represent the template fits derived by Hamuy et al. (1996c) to 
the corrected photometry. The value of E(B — V)Tempiate is defined as the difference between the 
template fits and the unreddened Lira relation (the dashed lines) at ty = 53 days (ts — 55 days). 

Fig. 4. — Comparison of color excesses calculated via the "Direct" and "Template" methods from 
the B — V evolution during the phase interval 30 < ty < 90. 

Fig. 5.— (Above) The B 

max ~ Vmax "color" is plotted versus the decline rate parameter Ami5(i?) 
for 20 events with E(B — V)Taii < 0.06. These SNe are likely to have little or no host galaxy 
reddening. (Below) The V max — I ma x "color" is plotted versus the decline rate parameter Amis(i?) 
for the same subsample of essentially unreddened SNe la. 

Fig. 6. — Comparison of color excesses calculated from the B max — V ma x and Vmax — Imax 
colors versus E(B — V)T a u- The solid line plotted in the lower panel assumes E(V — I)Max = 
1.25E(B — V)Max-, typical of dust in our own Galaxy. Open symbols correspond to SNe for which 
the first observation was not obtained until > 7 days after the epoch of B max . 

Fig. 7. — The total dust reddening (i.e., the sum of the Galactic and host galaxy contributions) 
plotted as a function of redshift for the 29 SNe la in the Calan/Tololo sample. The dashed line shows 
the upper limit to the allowed color excess for the limiting magnitude of this search, B Limit = 19. 

Fig. 8. — (Left) The absolute BVI magnitudes corrected only for Galactic reddening plotted versus 
the decline rate parameter Amis(i?) for the full sample of 41 Calan/Tololo and CfA SNe la with z 
> 0.01. (Right) The same diagram after elimination of the 23 SNe in the sample for which we find 



significant host galaxy reddening. The dashed lines correspond to the linear fits given by |Hamuy 



et al. (1996a) . 
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Fig. 9. — Hubble diag rams for the 26 SNe let in the Calan/Tololo "low-extinction" (£> max — Vmax ^ 
0.2) sample. The data are shown corrected for Galactic reddening only (upper panel), Galactic 
reddening and the Am 15 (i3) versus M max relation (middle panel), and Galactic reddening, host 
galaxy reddening, and the Ami5(i3) versus M max relation (lower panel). The Hubble constants 
and dispersions derived for each case are indicated. 

Fig. 10. — SN la Absolute magnitudes plotted versus B max — V max for two samples consisting of 
the Calan/Tololo + CFA SNe la with z > 0.01 and nearby SNe la for which distances have been 
determined via Cepheids. Distances to the former set of objects were calculated assuming H Q =63.3. 
The decline rate dependence of both the absolute magnitudes and the B max — V max color have been 
removed via equations 20-2^. For both samples, only those SNe with light curve coverage beginning 



< 5 days after the epoch of B max were considered. The canonical Galactic reddening vectors are 
illustrated as solid lines. 
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Table 1. Photometry References for Nearby SNe la 



SN 


References 


1937C 


Pierce & Jacoby (1995) 


1972E 


Ardeberg & de Grood (1973); Cousins (1972); Eggen & Phillips (unpublished); Lee et al. (1972) 


1980N 


Hamuy et al. (1991) 


1981B 


Barbon et al. (1982); Buta & Turner (1983); Tsvetkov (1982) 


1986G 


Phillips et al. (1987) 


1989B 


Wells et al. (1994) 


1990N 


Lira et al. (1998) 


1991T 


Lira et al. (1998) 


1991bg 


Filippenko et al. (1992); Leibundgut et al. (1993) 


1992A 


Suntzeff et al. (unpublished) 


1994D 


Smith et al. (unpublished) 


1996X 


Covarrubias et al. (unpublished) 


1998bu 


Suntzeff et al. (1999) 
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Table 2. Host Galaxy Reddenings 



DM 


A T-v-i ( T3\ a 
Am 15l^J fe s 


t,(±S - V )Gal 


TiV R T/"\ b 
*H-D - V frail 


TiV R T/"\ b 
h,{B - V )Max 


T7YT/ J\ b 
t,(V - l)Max 


— V )Avg 


1937C 


0.87(10) 


0.014(001) 


0.012(053) 


0.061(060) 




0.03(03) 


1972E 


0.87(10) 


0.056(006) 


-0.003(051) 


0.040(060) 


-0.001(073) 


0.01(03) 


1980N 


1.28(04) 


0.021(002) 


0.110(052) 


0.077(044) 


-0.009(056) 


0.05(02) 


1981B 


1.10(07) 


0.018(002) 


0.084(053) 


0.150(044) 




0.11(03) 


1986G 


1.73(07) 


0.115(011) 


0.668(057) 






0.50(05) 


1989B 


1.31(07) 


0.032(003) 


0.457(067) 


0.363(069) 


0.465(077) 


0.34(04) 


1990N 


1.07(05) 


0.026(003) 


0.160(057) 


0.106(051) 


0.057(061) 


0.09(03) 


1990O 


0.96(10) 


0.093(009) 


0.069(060) 


0.072(061) 


-0.071(073) 


0.02(03) 


1990T 


1.15(10) 


0.053(005) 


0.138(056) 


0.080(104) 


0.052(112) 


0.09(04) 


1990Y 


1.13(10) 


0.008(001) 


0.348(061) 


0.388(106) 


0.064(111) 


0.23(04) 


1990af 


1.56(05) 


0.035(004) 


-0.008(111) 


0.052(047) 




0.04(03) 


1991S 


1.04(10) 


0.026(003) 


0.122(069) 


0.089(104) 


0.004(111) 


0.06(04) 


1991T 


0.94(05) 


0.022(002) 


0.174(052) 






0.14(05) d 


1991U 


1.06(10) 


0.062(006) 


0.145(063) 


0.141(105) 


0.152(112) 


0.11(04) 


1991ag 


0.87(10) 


0.062(006) 


0.075(054) 


0.122(061) 


0.030(080) 


0.07(03) 


1991bg 


1.93(10) 


0.040(004) 


0.038(064) 






0.03(05) 


1992A 


1.47(05) 


0.017(002) 


0.004(051) 


0.000(046) 


-0.061(059) 


0.00(02) 


1992J 


1.56(10) 


0.057(006) 


0.050(065) 


0.129(106) 


-0.064(114) 


0.03(04) 


1992K 


1.93(10) 


0.101(010) 


-0.029(055) 






0.00(04) 


1992P 


0.87(10) 


0.021(002) 


0.082(058) 


0.064(046) 


0.093(072) 


0.07(03) 


1992ae 


1.28(10) 


0.036(004) 


0.215(097) 


0.132(060) 




0.12(04) 


1992ag 


1.19(10) 


0.097(010) 


0.622(066) 


0.151(061) 


0.116(072) 


0.10(04) c 


1992al 


1.11(05) 


0.034(003) 


0.052(053) 


-0.015(044) 


-0.013(068) 


0.01(03) 


1992aq 


1.46(10) 


0.012(001) 


0.137(096) 


0.115(061) 


-0.200(074) 


0.00(04) 


1992au 


1.49(10) 


0.017(002) 


0.031(093) 


0.058(105) 


-0.243(113) 


0.00(04) 


1992bc 


0.87(05) 


0.022(002) 


-0.055(053) 


-0.006(044) 


0.012(069) 


0.00(02) 


1992bg 


1.15(10) 


0.185(018) 


0.050(059) 


0.009(062) 


-0.016(075) 


0.01(03) 


1992bh 


1.05(10) 


0.022(002) 


0.177(071) 


0.132(060) 


0.126(071) 


0.12(03) 


1992bk 


1.57(10) 


0.015(001) 


0.034(065) 


0.001(062) 


-0.013(114) 


0.01(03) 


1992bl 


1.51(10) 


0.011(001) 


-0.006(065) 


0.012(061) 


-0.073(075) 


0.00(03) 


1992bo 


1.69(05) 


0.027(003) 


-0.004(058) 


-0.014(049) 


0.020(064) 


0.00(03) 


1992bp 


1.32(10) 


0.069(007) 


0.145(077) 


-0.034(060) 


-0.090(073) 


0.00(03) 


1992br 


1.69(10) 


0.026(003) 


-0.002(102) 


0.016(063) 




0.01(04) 


1992bs 


1.13(10) 


0.012(001) 


0.142(102) 


0.123(060) 




0.10(04) 


1993B 


1.04(10) 


0.079(008) 


0.200(069) 


0.196(061) 


0.047(072) 


0.12(03) 


1993H 


1.69(10) 


0.060(006) 


0.064(057) 






0.05(04) 
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Table 2 — Continued 



SN 


Am 15 (B) obs a 


E(B - V) Gal h 


E(B - V) Tail h 


E(B - V) M ax h 


E(V - I) M ax h 


E(B - V) A vg & 


19930 


1.22(05) 


0.053(005) 


0.059(059) 


-0.046(044) 


0.025(069) 


0.00(03) 


1993ac 


1.19(10) 


0.163(016) 


0.080(089) 


0.096(062) 


0.250(074) 


0.12(04) 


1993ae 


1.43(10) 


0.039(004) 


-0.024(053) 


-0.036(104) 


-0.009(113) 


0.00(03) 


1993ag 


1.32(10) 


0.112(011) 


0.152(061) 


0.112(061) 


-0.025(074) 


0.07(03) 


1993ah 


1.30(10) 


0.020(002) 


0.075(074) 


-0.013(104) 


-0.062(112) 


0.01(04) 


1994D 


1.32(05) 


0.022(002) 


-0.060(052) 


-0.017(044) 


-0.010(057) 


0.00(02) 


1994M 


1.44(10) 


0.024(002) 


0.115(064) 


0.057(061) 


0.126(074) 


0.08(03) 


1994Q 


1.03(10) 


0.017(002) 


0.090(065) 


0.099(104) 


0.048(111) 


0.06(04) 


1994S 


1.10(10) 


0.021(002) 


0.001(112) 


0.009(045) 


-0.034(071) 


0.00(03) 


1994T 


1.39(10) 


0.029(003) 


0.105(114) 


0.126(061) 


0.112(073) 


0.09(04) 


1994ae 


0.86(05) 


0.031(003) 


0.184(056) 


0.175(045) 


0.053(057) 


0.12(03) 


1995D 


0.99(05) 


0.058(006) 


0.103(054) 


0.044(044) 


0.016(056) 


0.04(02) 


1995E 


1.06(05) 


0.027(003) 


0.934(072) 


0.793(045) 


0.992(068) 


0.74(03) 


1995ac 


0.91(05) 


0.042(004) 


0.109(054) 






0.08(04) d 


1995ak 


1.26(10) 


0.043(004) 


0.260(057) 


0.107(060) 


0.294(072) 


0.18(03) 


1995al 


0.83(05) 


0.014(001) 


0.157(064) 


0.246(045) 




0.15(03) 


1995bd 


0.84(05) 


0.495(050) 


0.242(074) 






0.15(06) d 


1996C 


0.97(10) 


0.014(001) 


0.078(055) 


0.139(060) 


0.106(072) 


0.09(03) 


1996X 


1.25(05) 


0.069(007) 


-0.010(051) 


-0.016(044) 


0.076(056) 


0.01(02) 


1996Z 


1.22(10) 


0.063(006) 


0.500(073) 


0.352(061) 




0.33(04) 


1996ai 


0.99(10) 


0.014(001) 


2.015(080) 


1.936(065) 


1.575(072) 


1.44(04) 


1996bk 


1.75(10) 


0.018(002) 


0.274(064) 






0.19(05) 


1996bl 


1.17(10) 


0.105(011) 


0.113(063) 


0.095(061) 


0.090(073) 


0.08(03) 


1996bo 


1.25(05) 


0.078(008) 


0.358(053) 


0.385(061) 


0.253(069) 


0.28(03) 


1998bu 


1.01(05) 


0.025(003) 


0.328(053) 


0.375(054) 


0.519(074) 


0.33(03) 


1996bv 


0.93(10) 


0.105(010) 


0.215(062) 


0.233(061) 


0.318(073) 


0.21(03) 



a The mean errors are listed in units of CPOl. 
b The mean errors are listed in units of CP 001. 

C E(B — V)Taii measurement ignored in calculation of E(B — V)A V g due to liklihood of systematic errors 
in late-time photometry. 

d 1991T-like event; E(B - V) Avg based only on value of E(B - V) Ta il- 
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Table 3. Fits to Decline Rate versus Luminosity Relation a 



AM max 


= a[Am 15 (B) 


- 1.1] + 6[Ami5(B) - 


l.l] 2 




Bandpass 


a b 


b b a (mag) 


"X.V 


n 


B 


0.786(398) 


0.633(742) 0.11 


0.47 


17 


V 


0.672(396) 


0.633(742) 0.09 


0.42 


17 


I 


0.422(400) 


0.633(742) 0.13 


0.83 


15 


a AM max = 


M max M maa 


; (Ami 5 (B) = 1.1) 






b The errors 


are listed in units of 0™001. 
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E(B-V) Template 
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CD 
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-0.4 
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a 0.0 



l 



_ E 0.2 



0.4 



0.6 




• < 7 days after B max 
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- MLCS (Riessetal. 1996) 

- New MLCS (Riess et al. 1998) 
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Calan/Tololo "Low Extinction" Sample 
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